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ABSTRACT 

We explore here the young stellar populations in the Rosette Molecular Cloud 
(RMC) region with high spatial resolution X-ray images from the Chandra X-ray 
Observatory, which are effective in locating weak-lined T Tauri stars as well as 
disk-bearing young stars. A total of 395 X-ray point sources are detected, 299 of 
which (76%) have an optical or near-infrared (NIR) counterpart identified from 
deep FLAMINGOS images. From X-ray and mass sensitivity limits, we infer a 
total population of ~ 1700 young stars in the survey region. Based on smoothed 
stellar surface density maps, we investigate the spatial distribution of the X-ray 
sources and define three distinctive structures and substructures within them. 
Structures B and C are associated with previously known embedded IR clusters, 
while structure A is a new X-ray-identified unobscured cluster. A high mass 
protostar RMCX #89 = IRAS 06306+0437 and its associated sparse cluster 
is studied. The different subregions are not coeval but do not show a simple 
spatial-age pattern. Disk fractions vary between subregions and are generally 
< 20% of the total stellar population inferred from the X-ray survey. The data 
are consistent with speculations that triggered star formation around the HII 
region is present in the RMC, but do not support a simple sequential triggering 
process through the cloud interior. While a significant fraction of young stars are 
located in a distributed population throughout the RMC region, it is not clear 
they originated in clustered environments. 

Subject headings: ISM: clouds - ISM: individual (Rosette Nebula) - stars: forma- 
tion - stars: pre-main sequence - X-Rays: stars 



department of Astronomy & Astrophysics, The Pennsylvania State University, 525 Davey Lab, Univer- 
sity Park, PA 16802; jwang@astro.psu.edu, edf@astro.psu.edu, townsley@astro.psu.edu 

2 Centro Astronomico Hispano Aleman, Camino Bajo Huetor 50, Granada, Spain 18008 

3 Department of Astronomy, University of Florida, 211 Bryant Space Science Center, Gainesville, FL 
32611 



-2- 



Introduction 



In the last two decades, developments in infrared (IR) technologies, such as JHK arrays 



processes in molecular clouds (e.g., 


Lada et al. 


1991 




Lada k, Lada 


1995 




Carpenter 


2000; 


Evans et al. 


2003: 


Churchwell et al. 


2004; 


Young et al. 


boos 




Harvev et al. 


2006). One ma- 



jor discovery from such surve ys of nearby clouds i s that a large fraction (~70%-90%) of stars 
form in embedded clusters ( jLada fc Ladal 120031 ). The spatial distribution and hierarchi- 
cal structure of embedded clusters in giant molecular clouds (GMCs) should reveal physical 



processes in clu s ter formation (e.g., Elmegreen |2000l ; iBonnell et al.ll200ll ; lLada fc Lada 



Tan et al. 2006; Ballesteros-Paredes et al. 



2003 



20071 ). It is particularly important to understand 



the influence of initial clusters on later star formation in adjacent cloud material. The expan- 
sion of the HII region around massive stars in young stellar clusters can, i n principle, sequen- 
tially trigger for mation of second-generation OB clusters , as described by Elmeg reen fc Lada 



(119771 ) (see also IWhitworth et al.lll994l ; bale et all 120071 : IPreibisch fc Zinneckerl 120071 ). How- 



ever, the operation of triggering in specific cases of OB associations in GMCs is often unclear. 

As one of the most massive Galactic GMCs, the Rosette Molecular Cloud (RMC) is ideal 
for studying the f ormation and po ssible triggering of embedded st ar clusters. With a tota l 
extent of ~150 pc JCox et al.l |l99oh and ~ 10 5 M^ of molecu lar gas (blitz fc Thaddeusl ll98oh. 
it lies only 1 .4 kpc from the Sun (IHensberge et al.ll2000l ). iTownsley et al.l (120031 . hereafter 



TFM03) and lRoman-Zuniga et al. (l2008bl . hereafter RL08) review past and present research 
on this well-studied and important star formation region. Its orientation is perpendicular 
to the line-of-sight with a large, annular HII region at its tip known as t he Rosette Nebula . 
This HII region is powered by the ~ 2 Myr old OB cluster NGC 2244 JWang et all I2008I . 
hereafter Paper I). The molecular gas exhib i ts the strongly clu mped structure characteristic 
of turbulent GMCs Jwilliams fc Blitz! Il998l : iHever et alJbood ). 



Phelps fc Ladal (119971 . henceforth, PL97) identified seven embedded clusters in the 
dense molecular cloud cores from stellar concentrations in near-IR (NIR) JHK images 
which are associated with IRAS mid-IR sources (see Figure la in Paper I). They suggested 
that both triggered and spontaneous star formation may be contributing to the forming 
clusters in this cloud . The recent deeper FLAMINGOS NIR imaging survey of Rosette 
( IRoman-Zuniga et al.l l2008al . hereafter REFL08) reveal four more clusters across the com- 
plex. They find an increasing fraction of younger stars (identified from i^-band excesses 
associated with hot circumstellar disks) with increasing distance from NGC 2244, implying 
a temporal sequence of star formation across the complex. Four new clusters and over 700 
stars with dusty disks are located in mi d-IR maps obtained with the IRAC and MIPS de- 
tectors on the Spitzer Space Telescope ( jPoulton et al.l 120081 ) . Both the FLAMINGOS and 
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Spitzer studies argue that the embedded star formation is not principally produced by se- 
quential triggering from the NGC 2244 OB cluster. 

Lying 2° from the Galactic Plane, IR imaging surveys in the Rosette region are unavoid- 
ably contaminated by the overwhelming field star populations. Therefore, identification of 
young embedded populations has relied on selecting stars with IR color excesses produced 
by dusty circumstellar disks. Because the many diskfree pre-main sequence stars (PMS) 
(weak-lined T Tauri stars, WTTS) that are cluster members cannot be identified in this 
way, statistical subtraction of the contaminating field population is needed to estimate the 
full embedded population from IR images. Even in very young clusters, the WTTS may 



outnumber stars with prominent disks (e.g.. IWalter et al.lll988t iNeuhaeuser et al.lll995l ). 



PMS stars both with and without disks are well- known to emit X-rays that can penetrate 
the heavy extinction of GMCs (IFeigelson et al.ll2007l ). The level of X-ray emission from young 
stellar populations is elevate d 10 — 10 4 times abov e the levels of X-ray emission from the 
old Galactic disk population (IPreibisch et al.ll2005l ). As a consequence, X-ray observations, 
particularly with the high spatial resolution Chandra X-ray Observatory, are highly efficient 
in selecting PMS members in stellar clusters. Chandra images provide remarkably complete 
and unbiased PMS populations with subarcsecond positions. X-ray surveys thus complement 
IR surveys in establishing the census of embedded PMS populations in GMCs which in turn 
can address issues relating to triggered star formation. Chandra has made recent advances 
in the study of star formation triggering in the vicinity of HII regions produced by the Cep 
OB3b, Trumpler 37 (IC 1396 ), NGC 6618 (M 17) and Tmmpler 16 (Carina) OB associations 
and a high latitude cluster (IGetman et al.ll2006l 120071 ; iBroos et al.l 120071 ; ISanchawala et al. 
20071 : iGetman et all [2008^ 1 . 



A previous X-ray study of the RMC was presented by iGregorio-Hetem et al.l (119981 ) 
based on an image from the ROSAT satellite. They report faint X-ray sources associated 
with embedded IR-emitting low mass star clusters. However the poor spatial resolution and 
soft energy band of ROSAT prevents identification of individual embedded stars. The diffuse 
and stellar emission from the central Rosette Nebula cluster NGC 2244 were studied in detail 
in TFM03 and Paper I. 

We now examine a mosaic of four Chandra images in the RMC, presented without de- 
tailed analysis in TFM03, to identify the PMS population and study the cluster formation 
process in molecular clouds, complementing the NIR studies (PL97, RL08, REFL08). Fig- 
ure 1 outlines the location of our images with respect to the molecular gas and dust emission 
and the NIR clusters. Section [2] describes the observations with § 12.21 outlining the stellar 
counterparts of the X-ray sources. Important individual stars are discussed in § [3j Charac- 
terization of the embedded clusters, and non-clustered (distributed) population, follows in 
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§ HI In particular we discover a previously unidentified older cluster (§ 14.41) and investigate 
clustered star formation in the particular environment of the Rosette Complex (§ 14.81) . In 
§ we summarize our findings and discuss the implications for star formation modes in the 
RMC. 



2. X-ray Sources in the RMC 



2.1. Chandra Data Analysis 

The RMC was observed with the Imaging Array of Chandra's Advanced CCD Imaging 
Spectrometer (ACIS-I) in a sequence of 20 ks exposures ob tained in January 2001 (ObsIDs 
1875-77). The observatory and instrument are described by IWeisskopf et al.l (120021 ). Details 
on the data filtering, point source detection and extraction, variabilit y analysis, and spe ctral 
fitting procedures using the ACIS Extract (AE) software packag^ ( Broos et al. 2002 ) are 
given in Paper I. This paper also defines the quantities presented in the source tables below. 

The resulting RMC source list contains 395 X-ray detections which we designate RMCX 
1 through 395. We divide these into a primary list of 347 highly reliable sources (Tabled]) and 
a secondary list of 48 tentative sources (Table [2]) with Pb > 0.001 likelihood of being spurious 
background fluctuations based on Poisson statistics. This approach ha s been adopted in our 



previous papers on Chandra observa tions of star forming regions (e.g. iTownsley et al.l 12006 
Wang et al.l 120071 ; iBroos et al.l 120071 ) ; it allows the reader to evaluate the validity of faint 
sources individually. Tableland Table [2] have a format that is identical to Tables 1 and 2 in 
Paper I. A Chandra mosaic of the RMC region (0.5-8 keV) overlaid with source extraction 
regions (calculated with AE) is shown in Figure [2j The image here is shown at reduced 
resolution (2" per pixel). 

During the ~ 20 ks observations, eleven sources display significant variability (Prs < 
0.005, shown in column 15 of Tables[T]and [21 where Pks is the significance of the statistic for 
1-sided Kolmogorov-Smirnov test comparing a uniform count rate model to the distribution 
of source event time stamps). These are associated with magnetic reconn ection flares which 
are responsible for most of PMS X-ray emission (jFeigelson et al.ll2007l ). Only one of the 
variable sources (source #119) has more than 100 net counts due to the short exposure; its 
light curve shows a slow 3-fold rise in emission over 2 hours followed by a decay over the 
subsequent 3 hours. Such slow-rise X-ray flares are not infrequently seen from PMS stars 



1 The ACIS Extract software package and User's Guide are available online at 
http : //www. astro .psu. edu/xray/docs/TARA/aejusers_guide .html 
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f lGetman et al.ll2008bh . 



Spectral analysis results for 164 sources with photometric significance Signif > 2.0 
(Signif defined as the net counts value divided by the error on that value; see Table [T] foot- 
note) are presented in Table [3j The table entries are described in the table footnot es. ACIS 



spectral distributions were fit to thermal plasma models ("apec:".ISmith et al.ll2001l) with in 



terstellar soft X-ray absorption ("wabs;". iMorrison fc McCammonlll983l ) convolved with the 
telescope mirror and detector spectral response. Best-fit absorbing column densities range 
from negligible to l og Nh ~ 23.5 cm -2 , equivalent to a visual absorption of Ay ~ 200 mag 
(jVuong et al.ll2003l ). Temperatures range from kT ~ 0.3 keV to the hardest measurable by 
ACIS which we truncate at kT = 15 keV. The range of total band (0.5 — 8 keV) absorp- 
tion corrected luminosities (L tlC ) derived from spectral modeling is 29.6 < \ogL t c < 32.4 
ergs s _1 . Assuming a 2 keV plasma temperature and an average Ay = 5 mag visual extinc- 
tion (log Nh ~ 21.9 cm -2 ), the faintest on-axis detection in Table [2] corresponds to a limiting 
luminosity of \ogL tiC ~ 29.3 ergs s _1 . Note however that the patchy extinction throughout 
our fields will change the limiting sensitivity of the observation somewhat (§ 12.31) . The 



brightest source in the field is source #119 with 369 ACIS counts in the 20 ks exposure. 
This count rate is not high enough to cause photon pile-up in the detector. 



2.2. Identification of Stellar Counterparts 



The X-ray sources, which are generally located with precision better than 0.5" in the 
2MASS/Hipparcos reference frame, are associated with counterparts at other wavelengths 
to establish their stellar properties. Optical coverage of the RMC in the literature is poor 
because of the high visual extinction towards the molecular cloud and the contamination 
by foreground stars. Few spectral types are measured. We ther efore use optical a nd IR 
photometric catalogs including: Whole-Sky USNO-B1.0 Catalog (IMonet et al.ll2003l . here- 
after USNO); 2MASS All-Sky Catalog of Point Sources JCutri et alil2003l . 2MASS): and 
the University of Florida FLAMINGOS Survey of Giant Molecular Clouds (PI: E. Lada). 
The FLAMINGOS observations of the Rosette Complex fields and IR data reduction are 
described by Roman-Zuniga (2006) and REFL08. We have also associated our sources with 
the coordinates of >50,000 stars seen in Spitzer IRAC 3.6/xm images (Poulton et al. 2008, 
private communication). Positional coincidence criteria are used to associate ACIS X-ray 



sources with optical and NIR (ONIR) sources as described in the Appendix of iBroos et al. 
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Associations between ACIS sources and ONIR sources are reported in Table HI Of the 
395 ACIS sources in the RMC, 299 (76%) have an ONIR counterpart identified. JHK magni- 
tudes from FLAMINGOS photometry are reported if available; we adopt 2MASS photometry 
for six bright stars (RMCX #25, 65, 164, 185, 209, 342) that are saturated in FLAMINGOS 
(H < 11 mag, Roman-Zuniga 2006). NIR counterparts of 279 ACIS sources have measure- 
ments in all three JHK bands. The SIMBAD and VizieR catalog services were examined for 
complementary information; these results are reported in the table f ootnotes. Useful sou rces 
include the TASS Mark IV Photometric Su rvey of the Northern Sky (Droege et al.| 2006) and 
two large optical surveys of the Rosette (IMassey et al.l Il995l ; iBerghofer fc Christian! |2002j) 
that partly covered our Chandra fields. Source #65 is too bright (J = 9.83) for its spectral 
type (GO) at the RMC's distance, leading to its classification as a foreground star. The rest 
of the stars with known spectral types show colors and magnitudes consistent with being 
intermediate- to high-mass cloud members. 

We have examined the spatial distribution of the 96 X-ray sources that do not have 
associated ONI R counterparts in the context of Digital Sky Survey optical plates and CO 
emission maps (IHeyer et al.l 120061 ). Over 20 of these sources are located in the brightest 
nebulosity of the photodissociation region (PDR) of the Rosette Nebula, where source de- 
tection becomes extremely difficult at ONIR wavelengths. We suspect that most of these 
are cloud members. About 50 others are concentrated in regions where the dense molecular 
cores are located and show X-rays with harder median energy (MedE). Most of these are 
probably heavily obscured new members of the embedded clusters, although a few may be 
extragalactic sources coincidentally lying behind the embedded clusters. 

Only a small fraction of the X-ray sources, perhaps 10 of the 96 that do not have ONIR 
counterparts, seem to be randomly distributed across the fields with a surface density consis- 
tent with the extragalactic background of active galactic nu clei. Extragalactic contaminants 



are seen in all Chandra fields of young stellar clusters (e.g.. lGetman et al.ll2006l ; IWang et al. 
20071 ; iBroos et al.ll2007l ). Considering the contamination analysis of these earlier studies, the 
absorption by the RMC gas, and the short 20 ks exposures considered here, we estimate a 
contamination of ~3% from extragalactic sources and ~ 1% from older Galactic field star 
populations among the 395 ACIS sources, and even fewer among the sources with ONIR 
counterparts listed in Table HI 



2 Software implementing the matching algorithm is available in the TARA package at 
http : //www. astro .psu. edu/xray/docs/TARA/ 
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2.3. Limiting sensitivity and completeness 



The sensitivity limits of our X-ray observations are spatially uneven due to the patchy 
absorption in the RMC and to degraded point spread function (PSF) of the Chandra mirrors 
at large off-axis angles. Ba sed on previous experience with wavelet source detection and 
the ACIS Extract software iGetman et all bood : IvVang et all l2007i : iBroos et al1l2007h . we 
estimate that the far-off-axis sensitivity is approximately 2 times worse than on-axis (see 
§ 14.31) . Sources with 8 or more net counts would be detected anywhere in the three fields; 
over 60% of the sources in Table [1] lie in this complete regime. 

The typical extinction is Ay ~ 5 — 10 in the RMC regions with molecular concentrations 
and Ay ~ 1 in less obscured regions (REF L08). Assuming a kT = 2 keV thermal plasma 
as the spectral model of a typical PMS star (IPreibisch et al.ll2005l ). the absorption-corrected 
luminosity limit of 8 ACIS counts in 20 ks in the hard band (2-8 keV) is log Lh, c ~ 29.6 erg 
s^ 1 in obscured regions and ~ 30.0 erg s _1 in less obscured regions of the RMC (the observed 
luminosity in the 0.5-8 keV band logL t ~ 29.7 erg s -1 and ~ 30.1 erg s" 1 , respectively) 



Due to the strong empirical correlation between L x and PMS stellar mass (ITelleschi et al. 
20071 ). these luminosity limits correspond to a mass limit around ~ 0.5M Q . 



Adopting the X-ray luminosity function (XLF) of the Chandra Orion Ultradeep Project 
(COUP) observation of the Orion Nebula Cluster as a calibrator for other PMS stellar 
populations, we can estimate the completeness of o ur RMC observa t ions w ith respect to the 
full underlying PMS population. Using the XLF of iFeigelson et al.l (120051 ). we find that the 
RMC observations probe ~ 20% (~ 30%) of the PMS population down to the stellar limit 
of M = 0.08 M in the obscured (lightly obscured) regions of the RMC. These fractions are 
adopted to obtain estimates of the total stellar populations in various RMC structures below 
(§land Table E]). 



3. X-rays from Previously Studied Stars 



3.1. The 07 Star HD 46485 



The only previously known O star in our fields, HD 46485 (=BD+04°1318, 07V; 
WalbonJE97a Mafz-Apellaniz et al.ll2004l ) is detected with 310 counts in the 20 ks exposure 
(RMCX #164). A single-temperature thermal plasma model gives an acceptable fit to its 
X-ray spectrum with kT ~ 0.7 keV and a low absorbing column Nh ~ 2.1 x 10 21 cm -2 . 
The fit can be improved with a two-temperature thermal plasma model (Figure HJ), with 
a soft component {kT\ ~ 0.3 keV), a second harder component (kT 2 ~ 0.9 keV), and an 
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absorbing column of Nh ~ 4.1 x 10 21 cm -2 {Ay ~ 2.5). The soft plasma and the derived 
X-ray luminosity, logL t)C ~ 32.1 erg s -1 , are similar to the NGC 2244 O stars (e.g., HD 
46223, HD 46202; Paper I), consistent with being associated with instabilities in the radia- 
tively accelerated wind. Note that some O stars exhibit a much harder X-ray component; 
for e xample, 9 1 Ori C ( 07V) ionizing the Orion Nebula has a component with kT ~ 3 keV 
(e.g. 



Gagne et al.ll2005l ). This hard component is attributed to a magnetically channeled 



wind shock (MCWS). These differences in the X-ray properties of O stars of similar spectral 
types may be attributed to the various magnetic field strength in different stars; in this 
case, HD 46485 lik ely does not have a strong field. Mea surements of magnetic properties of 



massive stars (e.g.. lAlecian et al.ll2008t iPetit et al.ll2008l ) will be able to further evaluate the 



link between X-ray emission and the presence of magnetic fields in OB stars. 

HD 46485 does not appear to have any relationship to the embedded clusters of the 
RMC. It is much less absorbed and lies outside all of the clusters discussed here or in the 
infrared literature. The reported proper motio n of HD 46485 is smal l (/i Q = -0.26 mas/yr, fj,$ 
= -1.55 mas/yr, The HIPPA RCOS Catalogue; iPerryman et al.lll997f ). and the radial velocity 
is ~ 15 km s _1 (lEvanslll967l ). It might be a runaway O star from NGC 2244, or a somewhat 
older member of the open clusters in the larger OB association (jPerezlll99ll ). 



3.2. RMCX #89 = IRAS 06306+0437: A new massive protostar 



This star in PL97's cluster #2 (PL2) exhibits unusually red colors in the NIR, with 
J — H = 3.22 and H — K = 2.41. These colors imply both unusually strong reddening 
(Ay ~ 30) and a strong .fT-band excess. It is also a bright MSX point source (MSX6C 
G206. 7804-01. 9395). The Chandra source has only 21 net counts but the X-ray properties 
are all consistent with a highly absorbed source: the median photon energy i s very hard (4.6 



keV) which, based on the Orion Nebula calibration between MedE and Nh (IFeigelson et al. 
20051 ). corresponds to absorption log Nh = 23.2 cm -2 or Ay ~ 100. The observed X-ray 
luminosity in the hard band (2 — 8 keV) is logL^ ~ 31.2 erg s _1 . While not impossibly 
high for a lower mass PMS s tar, this is at the top ~ 1% of the PMS X-ray luminosity 



function (IFeigelson et al.l 120051 ) . Among protostars, only one of several dozen Class I in the 
Orion Nebula region flares up to a compa rable X-ray level (e.g., COUP #554 in OMC IS, 
Grosso et al.l 120051 ; iPrisinzano et al.l 120081 ) . 



Although the evolutionary stage of PMS stars cannot be directly established from hard- 
band X-ray properties because magnetic reconnection flares are present throughout the Class 
I-III phases, our past experience suggests that sources with absorption exceeding log Nh ~ 
22.3 cm -2 are protostars where most of the absorption arises in local envelopes rather than 
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the larger-scale molecular cloud (IGetman et al.l 120061 ). A protostellar classification is also 
indicated by the rapidly rising flux (high 24/8 /zm flux ratio) and roughly fla t 25 — 100 Aim 



spect ral energy distribution measured with the Spitzer and IRAS satellites ( IPoulton et al. 
2008h . 



RMCX #89 = IRAS 06306+0437 is also a massive system. With a dereddened K 
magnitude around K ~ 7, it has a bolometric luminosity associated with an early B or late 
O star (see the NIR color magnitude diagram in Paper I). As O stars generally have softer 
X-ray spectra than PMS stars, the correction from hard band to the total 0.5 — 8 keV band is 
larger; for the 06-07 star 9 l C Ori in the Orion Nebula Cluster, this correction is a factor of 
~ 6. If we assume the same correction, the estimated absorption-corrected total band X-ray 
luminosity is logL t c ~ 32.0 erg s _1 for RMCX #89. This source is similar to a candidate 
massive YSO found in the Chandra observation of the Pismis 24 cluster in the NGC 6357 



star for ming complex (|Wang et al.ll2OO70 . In fact, IRAS 06306+0437 is a known water maser 
source (IBrand et al.lll994l ). It was included in a CS(2-1) survey of I RAS sources with far-IR 
color characteris tics of ultra compact HII regions, but not detected ( iBronfman et al.lll996l ). 
Wu et al.l (120061 ) observed IRAS 06306+0437 among other water masers not associated with 
known HII regions or low mass YSOs and detected NH 3 (1,1) emission coincident with the 
IRAS /MSX infrared peak. They suggested the ammonia core may be a high mass protostellar 
candidate. 

The IRAS catalog notes that the MIR source is extended. We examine the dusty envi- 
ronment of IRAS 06306+0437 in Figure which we extracted from the archived IRAC and 
MIPS data from the Spitzer satellite. The star lies within a bright region of heated dust 
surrounded by an elongated complex of infrared filaments extending over several arcmin- 
utes which is probably associated with the larger photodissociation region of the Rosette 
Nebula. A smaller-s cale view of the dust emission around IRAS 06306+0437 is shown by 
Poulton et al.l (120081 ) . They note that this star does not have an accompanying concentration 
of young Class I or II stellar companions; i.e. that PL 2 is not a true stellar cluster. Figure [5] 
shows 25 X-ray selected stars within 1 pc of RMCX #89. We have classified these as Class I, 
II or III based on a preliminary examination of their NIR colors, including two other Class I 
protostars, three Class II systems, and ~ 20 Class III systems. This is consistent with the 
concentration of NIR i^-band excess stars found in PL 2 (Figure 7 in REFL08). Another 
candidate protostar, RMCX #72 (K = 10.98), is close to RMCX #89, lying ~ 40" to the 
east in a secondary dust condensation. This further indicates the youth of this region. 

We thus conclude that RMCX #89 is a massive proto-0 star lying in a sparse cluster 
of PMS stars, most of which are older Class III systems. The implied age offset — that the 
central O star is younger than the cluster of lower-mass PMS stars — is similar to that re- 



-10- 



cently inferred from Chandra study of the embedded W3 Main cluster (IFeigelson fc Townsley 
2008h . 



3.3. Intermediate-mass X-ray Stars 



Block et al.l (119931 ) list the five brightest .fT-band sources in the PL4 region. The two 
bright stars BG/IRS1 and BG/IRS2 (acronyms used in Block et al. 1993 = B G) form a close 



pair whose spectra show Br7 and are likely Herbig Ae/Be (HAeBe) stars (IHanson et al. 



1993|). RMCX #209 and #207 are matched to BG/IRS2 and BG/IRS4, respectively. Both 
have ~50 counts and the spectral fits yield kT ~ 3 keV with absorption log Nh ~ 22.2 
cm" 2 for RMCX #207 and logA^ ~ 22.5 cm" 2 for RMCX #2 09. This absorption column 
(Ay ~ 15 — 20) agrees with the Ak ~ 1.5 mag derived by IHanson et al.1 (119931 ). The 
absorption corrected luminosities are logL t c ~ 31.5 erg s _1 . Their hard X-ray spectra are 
consis tent with the high X-r a y temperatures (^ 20 MK) seen in X-ray surveys of HAeBe 
stars jHamaeuchi et al.ll2005l : Istelzer et al.l l2006h . RMCX #210 coincides with BG/IRS 5 
but has too few counts for spectral fitting. 

The brightest X-ray source in our Chandra fields is RMCX #119 with 369 ACIS counts. 
The X-ray spectrum can be roughly fit with a standard two-temperature thermal plasma 
model but there seems to be some abundance anomaly around 1 keV. The derived X-ray 
luminosity is logL t c ~ 31.9 erg s _1 , too high for most lower-mass PMS stars but consistent 
with the range of HAeBe systems. The stellar counterpart is bright, with V = 13.9, and 
was found not to have strong Ha emission in a follow-up study of ROSAT sources in the 
Rosette region faerghofer fc Christian! hopj ). Its location in the NIR color-magnitude dia- 
gram indicates an intermediate mass st ar, although possib ly with weak or absent accretion 



as suggested by the weak Ha emission (jCalvet et al 



2000). 



3.4. Protostellar sources of Herbig-Haro outflows 



Five features of excited gas seen in the optical band [SB] survey by lYbarra fc Phelps 



( 12004J ) lie in our Chandra field. These are probably Herbig-Haro outflows from protostars. 
We have carefully examined the X-ray fields in the vicinity of these sources. In a region 
6' southwest of PL 2, the X-ray emitting star RMCX #51 is found near the extended [SII] 
feature RMC-F. Its IR counterpart is a heavily reddened star with K = 12.94, which has no 
significant .fT-band excess. This star may be responsible for the [SII] emission feature. 



RMC-H is an arc-shaped feature of excited gas outlining a globule of gas northwest of 
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PL 2 with a possible HH outflow originating in the globule. Three X-ray sources are seen 
coincidently tracing the edge of the globule. One of them is the Class I protostar RMCX #58 
which may be powering the jet-like outflow. 



Ybarra fc Phelpd (12004 ) suggests that RMC-L is an HH outflow driven by a star in PL 4 
associated with IRAS 06314+0427. In this context, the X-ray detected Class I protostar 
RMCX #315 may be the exciting source as it lies eastwards on the axis of symmetry of 
the bow-shock feature. The 3.5 pc separation between the star and shock is rather large, 
but comparable large-scale collimated flows have been detected in other star forming regions 
(e.g., HH80/81, > 5 pc long; Marti et al. 1993). Other [SII] features such as RMC-E, 
RMC-G, and RMC-M do not have an obvious correlation with our X-ray sources. 



4. Stellar Clusters in the Chandra Population 

Visual examination of the spatial distribution of X-ray sources in Figure [2] shows a 
non-random clustering. This is also evident in the earlier presentation of the field shown in 
Figure 6a of TFM03. We identify the X-ray source concentrations, which can be considered 
to be 'star clusters,' using two methods in §4.1114.21 The newly identified XA cluster is 
describe in §4.41 and the previously established embedded XB and XC clusters are discussed 
in §4.51 Tables [Mni provide quantitative results. 



4.1. Clusters from the Smoothed X-ray Source Distribution 



Wang et al.l (120071 ) and iBroos et al.l (120071 ) found that smoothed maps of the stellar 



distribution derived from Chandra source lists are effective in revealing the structure of 
young stellar clusters. For example, in M 17, the X-ray stellar density maps highlighted the 
region's complex structure, including the compact NGC 6618 cluster, a triggered population 
along the edge of the HII region, several small embedded clusters, and a previously unknown 
subcluster. A similar approach applied to 2MASS star distrib utions revealed the NGC 2237 



subcluster on the western side of the Rosette Nebula (jLill2005l ). Compared to infrared maps, 
X-ray maps are less affected by contamination from Galactic field stars and show richer 
populations than disky IR-excess samples. 



Smoothed maps of the Chandra sources in the RMC are shown in Figure Ian. The three 



3 Note that each star has the same weight in these maps, so they do not accurately reflect the distribution 
of X-ray emission. A smoothed map of the X-ray emission is shown by TFM03. 
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ACIS fields have closely matched 20 ks exposures, only a small sensitivity bias is present 
between the on-axis and off-axis areas of each field, and contamination by extraneous sources 
is small ( §2.2p . The three rows of panels in Figure M show the total source sample, lightly 
obscured sources with MedE < 2.0 keV, and heavily obscured sources with MedE > 2.0 
keV. The left and right panels show maps with 3' and 2' radius (1.2 and 0.8 pc) smoothing 
kernels, respectively. 

Three large-scale structures are seen as overdensities in the total sample, irrespective 
of the adopted smoothing kernel size. We define these regions as XA, XB, and XC and 
trace their outline from a contour of the 3' kernel. We note seven small-scale overdensities 
within the three major structures in the 2' map of the total population. We label these 
XA1-XA3, XB1-XB3, and XC1. In addition, a distinct substructure of heavily obscured 
stars extending southward from XC1 is seen in panel Figure E]/ which we call XC2. Note 
that the substructures are defined from smoothed contours on different maps and are not 
contiguous; the star totals in the main structures are thus larger than the sum of the stars 
in the substructures. We label stars lying outside all overdensities "distributed" stars in the 
RMC. 

While the regions are defined based on subjectively chosen surface density contours, we 
are confident that these structures are not artifacts caused by the slightly degraded X-ray 
sensitivities at large off-axis angles on the ACIS-I CCDs. The overdensities are not located 
at the centers of the fields (where the telescope response is maximized) and represent a factor 
of 2 — 20 difference in stellar density which cannot be solely attributed to the instrumental 
sensitivity gradients. 



4.2. Clusters from the Nearest Neighbors Method 



A seco nd technique for identifying stellar cl usters is analysis of the distribution of nearest 
neighbors (jCasertano &: Hutlll985l ; lDigglell2003l ). This was used by REFL08 to establish the 
content of the seven PL97 clusters in the JHK catalog of the RMC. We apply the same 
nearest neighbor method used by REFL08 based on the local density derived from the 
distances to the 10 nearest neighbors to each star. 

Contours of the resulting nearest neighbor maps for our Chandra RMC observations are 
shown in Figure [7] for the total, unobscured, and obscured populations. The locations and 
morphologies of the resulting concentrations closely resemble the structures that we defined 
using the smoothed surface density maps for X-ray sources. This indicates the robustness of 
the inferred clustering irrespective of the methodology used to smooth the star distribution. 
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4.3. Evaluating the Sensitivity Across the Mosaic 

The maps shown in Figure [6] used all detected X-ray sources. However, the sensitivity 
across the field is not uniform; there is a bias towards detecting faint sources in the area 
close to the individual field aim points and in the overlapping regions with deeper exposure. 
It is possible that this sensitivity variations can create spurious structure in an intrinsically 
smooth distribution of sources. We evaluate the importance of this effect by creating a 
smoothed source map similar to Figure [6] which removes the sensitivity degradation due to 
telescope vignetting, effective exposure and background across the fields. 



Following iLehmer et al.l (120051 ) and iLuo et al.l (120081 ) , we estimate the sensitivity across 



the X-ray imaging mosaic assuming a Poisson model. The resulting relation can be approx- 
imately represented by logiV Km = 0.558 + 0.322(log6) + 0.13(log6) 2 + 0.037(log6) 3 , where 
Nu m is the required number of counts for a source required to be detected, and b is the 
background counts measured in a local region. The coefficients are obtained by constructing 
a lower envelope in the plot of Net Full counts against Bkgd Full counts tabulated in Table [TJ 
We then construct a sensitivity map across the mosaic by dividing a smoothed map of Ni im 
by the telescope exposure maps in units of cm 2 s (shown in Figure [H see also Appendix 
of TFM03). The resulting sensitivity map is typically 6 x 10 -7 counts cm -2 s _1 near the 
aimpoint and 1 x 10 -6 near the edge of a typical 20 ks field. It reaches 2 x 10 -7 in a small 
region around sub cluster XAl where the field overlaps the much deeper ~100 ks exposure 
of NGC 2244 discussed in Paper I. 

This corresponds to a ~8 counts limit across the 20 ks exposures and ~40 counts in the 
overlapping region around XAl. Applying this spatially varying counts threshold removes the 
sensitivity bias from telescope degradation and overlapping fields, and provides a complete 
sub-sample of our X-ray detections. This procedure removes 30% of the weaker sources in 
Table 1 (the remaining sources with number of counts above the detection threshold are 
marked in Table 1, see table footnote). The smoothed stellar surface density map of this 
reduced, spatially uniform sample is shown in Figure M 

The three major X-ray clusters, XA, XB, and XC, are present as prominent overden- 
sities. Some substructures are less significant than in the full sample map; XA2 and XB1 
are weakest with densities ~ 2 times above the background level. Additional possible sub- 
structures emerge: an enhancement at (6 /l 33 m 02 s , 4°47') east of XA3 and at (6 h 34 m 05 s , 
4°33') northwest of XC1. The reality of these features must be confirmed with deeper X-ray 
observations. 



-14- 



4.4. Identification of the new XA cluster 

The three main stellar enhancements and seven secondary enhancements are listed in 
Table [5] with approximate sizes. Most are around 3' or 1.2 pc in extent. Four of the seven 
substructures — XB2, XB3, XC1, and XC2 — are associated with clusters identified from 
JHK and Spitzer images (PL97, RL08, Poulton et al. 2008 ). Figure [T2l shows the outline 



of the three main structures on the molecular CO map of iHeyer et al.l (120061 ). It shows 
that these same four substructures-XB2, XB3, XC1, and XC2-are associated with dense 
molecular regions. 

X-ray cluster XA, on the other hand, lies in the evacuated region on the eastern side of 
the Ha-bright Rosette PDR. Located in the north-west side of the north-westernmost ACIS 
field where molecular material has been removed by the expanding HII region, its stars 
are less obscured compared to the XB and XC concentrations embedded in the molecular 
cloud. This is shown in columns 2 — 4 of Table El 30% of XA stars have MedE > 2 keV, 
corresponding to absorption \ogN H > 22.0 cm -2 or Ay > 6, compared to 41 — 46% in 
concentrations XB and XC. 

Most XA stars are concentrated in subcluster XA3 centered at (a, 5) = (06:32:38,-1-04:46: 1£ 
XAl and XA2 represent small extensions to the northwest and southwest, respectively. Sixty- 
one stars are in the XA X-ray sample. Scaling by the incompleteness factor for lightly ob- 
scured regions given in § 12.31 the estimated total stellar population in XA is 300 stars. This 
is about 15% of the total population estimated for the central NGC 2244 cluster using the 
same XLF-based method (~2000; Paper I). 

The lightly absorbed XA cluster was not clearly identified in previous IR surveys of the 
cloud. The location of substructure XAl is separated by ~ 20" from IRAS 06297+0453. 
However a close inspection of IRAS, DSS, and 2MASS images led us to conclude that this 
is a spatial coincidence: IRAS 06297+0453 is associated with the very bright foreground 
K0V star BD +04 1304 (V = 5.8, K = 3.69) and does not belong to the RMC. The XA 
cluster was not covered in the earlier PL97 survey. Six stars with .fT-band excesses are 
present in the FLAMINGOS survey (Table [6]), insufficient to highlight a cluster using the 
lOth-nearest neighbor method. However, it does appear as an enhancement in the smoothed 



surface density map of all 2MASS sources shown by|LjJ (120051 . Figures la and 2a), appearing 
as a stellar enhancement between XB and the main NGC 2244 cluster. A group of 32 stars 
with disky MIR spectral energy distributions is found within our XA c ontour in the Spitzer 



IRAC/MIPS, but again these were not identified as a distinct cluster by lPoulton et al.l (12008 
see their Figures 6b). Rather they appear as an eastward extension of the central NGC 2244 
cluster. The cluster is more spatially distinct from the main cluster in the ACIS X-ray mosaic 
of ~20ks pointings towards NGC 2244 and RMC (see Figure 6a in TFM03), which is not 
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sensitive to obscuration. 

The reason that the XA cluster emerges more clearly in our X-ray survey than in infrared 
surveys is probably due to presence of nebular contamination. It lies on a bright ridge of 



enha nced H Brackett-a and spatially variable PAH emission (Figures 1 and 4 of lPoulton et al. 



20081 ). These contaminants limit the detection or photometric accuracy of stars in certain 



infrared bands which then inhibits identification of disk emission from multiband spectral 
energy distributions. 



4.5. The Embedded Clusters XB and XC 



As noted above and in Table [5j the XB and XC concentrations of X-ray stars shown in 
Figures[6]and [7]correspond to embedded clusters established from IRAS, 2MASS, FLAMIN- 
GOS, and Spitzer studies. Their stellar populations inferred from NIR studies have been 
most thoroughly described by RL08. The XB and XC clusters have higher fractions of ob- 
scured stars than the XA cluster (46% and 41% vs. 30%; Table [6]). This, combined with 
their spatial relationship to dense molecular structures (Figure [T2l confirms that they are at 
least partly embedded in cloud material. 

The main advantage of our X-ray survey is that disk-free (PMS Class III) members are 
readily detected, subject only to the sensitivity limitations given in § 12.31 In most infrared 
surveys, Class III systems cannot be readily distinguished from the strong contamination of 
Galactic field stars. Infrared surveys are thus often limited to counting disky (PMS Class I 
and II) systems. The X-ray samples are thus usually larger than the IR-excess samples, and 
provide a less biased way to estimate disk fractions which is relevant to stellar ages. 

The XB cluster is the sparsest of the three X-ray concentrations with 52 detected stars 
and an estimated total population of 200 stars (Table [6]). It has an "L" -shaped structure with 
three s ubcomponents of sim ilar (~ 1 pc) sizes. Subcluster XB2 associated with Spitzer clus- 
ter C (IPoulton et al.l 120081 ) appears more embedded than the other subclusters. Subcluster 
XB3 is associated with NIR cluster PL 2 and MIR source IRAS 06306+0437. As discussed 
in § 13.21 above, this structure is dominated by the embedded massive star, RMCX #89. Both 
XB1 and XB2 contain bright MSX mid-IR sources (see footnotes to Table H]). 

The XC cluster is the richest of the three main X-ray structures in the cloud with 160 
detected stars. In the principal concentration of stars, XC1, roughly 2/3 are lightly obscured 
and 1/3 are heavily obscured. This rich clust er was first noted in if-band imaging by 
Block et al.l (119931). studied s pectro scopi callv in lHanson et al. Jl993f ). and noted NIR 
cluster by PL97, iBica et al.l (120031 ) , and iLi &: Smith! (120051 ) . The sparser XC2 component 
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is only evident in the stellar surface density maps for the obscured population (Figure [6^ 
and /) indicating it is deeply embedded. This region also appears much darker in the DSS 
optical image. It appears as a curved southern extension of the rich XC1 cluster in the 
X-ray maps, and is associated with the NIR clusters PL 5, REFL 8, and the Spitzer cluster 
E (Table [5]). The stellar structure may extend beyond the southeastern edge of our ACIS 
field. In contrast to XC1, only 1/4 of the stars in XC2 are lightly obscured as measured by 
the MedE indicator of soft X-ray absorption. XC2 is thus the most deeply embedded stellar 
cluster in the RMC regions we examine here. 



4.6. Evolutionary Stages of the Principal Clusters 

X-ray properties cannot discriminate very young clusters, dominated by accreting Class 
I-II systems, from older clusters, dominated by non-accreting Class III systems. The reason 
for this is well-established: the X-ray luminosity function and flaring properti es do not greatly 



change as low mass stars evolve from their Cl ass I through III phases (e.g. iPreibisch et al. 



20051 ; iTelleschi et all 120071 ; IStelzer et al.1 120071 ). X-ray selected samples are thus relatively 



unbiased with respect to age, and their ages can then be estimated from their rapidly evolving 
infrared properties. Here we focus on the high-quality JHK photometry available from the 
FLAMINGOS survey (REFL08) for X-ray sources within the boundaries of clusters XA, XB, 
and XC presented above. 

Figure [TU] shows the J — H vs. H — K color-color diagrams for the three clusters. 
The region between the left two dashed lines is associated with Class III objects (diskless 
WTTS) reddened by interstellar extinction. To the right of this band are sources exhibiting 
significant .fT-band excess; here we require E(H — K) > 2a(H — K) based on FLAMINGOS 
photometric errors for the sources to be K-band excess stars. The region between the middle 
and right-most dashed lines is occupied by Class II objects (PMS stars with circumstellar 
accretion disks), and protostars still possessing thick envelopes (Class I objects) lie beyond 
the right-most dashed line. Star counts in various regions are summarized in the middle 
columns of Table [6l 

The colors of most sources in cluster XA are consistent with PMS stars with small 
reddening, Ay ~ 1, typical of stars in the central NGC 2244 cluster illuminating the Rosette 
Nebula (Paper I). A few Class II sources are present showing higher Ay; these may lie within 
XA or behind XA in the molecular cloud. 



Sources in cluster XB show a wide range of absorptions; some sources are reddened 
with Ay > 5 and the .ff -excess sources are among the most reddened. This may suggest 
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the XB subclusters are distributed along the line of sight, at different cloud depths. Three 
Class I sources are found in this region. Despite the smaller number of sources in region XB, 
it has more IR-excess sources than the more populous XA region. Subcluster XB1 has 11 
detected X-ray stars with an inferred total population around 40 stars, but has no IR-excess 
systems. XB2 has more MIR-excess than NIR-excess stars, while XB3 has more NIR-excess 
than MIR-excess stars. We thus tentatively infer that the XB region is not homogeneous 
and coeval. Subcluster XB3, which includes RMCX #89 = IRAS 06306+0437 discussed in 
§ 13.21 is clearly the youngest structure in the region as the presence of NIR excess in the 
stars indicates inner disks have not been cleared yet. 

X-ray sources in cluster XC are mostly lightly obscured PMS stars with Ay ~ 2 and a 
dozen heavily obscured sources with Ay > 10. The XC cluster is very populous in the X-ray 
map, and only 5-10% of these stars have infrared excesses in the NIR or MIR bands. This, 
together with the presence of only one very young system in Figure [121 indicates that the 
region is not dominated by active star formation today despite the presence of two IRAS 
sources resolved into several IR clusters (PL 4, PL 5, REFL 8 and Spitzer's E). The very 
faint XC source RMCX #315 is associated with a faint NIR source with unusual colors, 
J — K > 3.8; it is either a Class I system or a very heavily reddened Class II system. NIR 
J vs. J — H color-magnitude diagrams shown in Figure [TT1 provide an approximate range of 
the mass distribution and local absorptions for the same stars shown in Figure [TUJ 

No candidate new OB stars appear in the RMC NIR color-magnitude diagram. This 
stands in marked contrast to Chandra studies of other prominent star forming regions where 
dozens of new obscure d OB stars were ide ntified by virtue of their X-ray em ission. These 
i nclude Westerlun d 1 Jskinner et al.ll2006h . and RCW 38 JWolk et al.l l2006h . NGC 6357 



flWang et all 120071 ). RCW 108 flWolk et al. 



20081 ) . There is no reason to believe that mas- 



sive star formation is suppressed in the RMC; more likely the population of the individual 
embedded cluster is simply too small to contain massive stars as expected from the initial 
mass function. 

We conclude that cluster XA, showing less absorption and very few IR-excess stars, 
appears older than clusters XB and XC. The latter structures are not unified, homogeneous 
clusters but rather collections of smaller stellar concentrations. Similar to the results of 
REFL08, we find no clear spatial pattern of cluster ages in this region. XC is the most popu- 
lous and, with only 5 — 10% fraction of IR-excess sources and no more than one protostar, is 
not very young. Active star formation is strongest in subcluster XB2 with several protostars 
including the massive system RMCX #89 = IRAS 06306+0437. 
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4.7. Relationship to the Molecular Gas Distribution 



Figure [12] shows the X-ray selected RMC stars, classified by evolutionary stage from 
the NIR color-colo r diagram (§ 14.61) . superposed on contours representing 12 CO emission 
( IHeyer et al.ll2006l ). The younger Class I and II sources are mostly seen around the CO 
molecular ridge near cluster XB and the molecular material near cluster XC. Class III stars 
appear more dispersed throughout the cloud and dominate the XA region where molecular 
material is no longer present. We can in terpret this tren d as reflecting the drifting of older 
Class III stars away from their birth sites (jFeigelsonll 1 9961 ) . For example, with a characteristic 
transverse speed of ~1 km s _1 , a star would travel 1 pc (~ 2.3') in 1 Myr. This is sufficient 
to disperse such stars originally formed in the clusters throughout the Chandra fields in a 
few Myrs, resulting in a "dispersed" population. 

The relationship of cluster XA to molecular gas is clearly different than for clusters XB 
and XC. XA lies within the bright nebulosity of the Rosette HII region where most of the 
molecular materials are gone. This cluster of less embedded and older PMS stars (see § 14.61) 
is likely an early population triggered by the expansion of the HII region. Its similarity to 
NGC 2237, a similar lightly obscured cluster on the western side of the Rosette Nebula, 
will be discussed in our Paper III (Wang et al., in preparation). The overdensity of stars in 
XB, including the NIR PL 2 cluster, is situated on the 'rim' of molecular material adjacent 
to the expanding Rosette PDR. It includes XB2, the youngest subcluster in our examined 
portion of the RMC. This cluster, and possibly the similar PL 1 cluster located south-west 
of our Chandra fields in a molecular clump, are reasonably interpreted as triggered by the 
expanding shell of the NGC 2244 HII region. The molecular clumps are probably undergoing 
photoevaporation and will eventually become unobscured clusters like region XA. 

Region XC contains the large and most extended population of embedded stars, includ- 
ing several infrared clusters. Its relationship to a dense molecular clump has been noted by 
Heyer et al.l (120061 ) and REFL08. Our analysis above (§ 14.61) suggests that subcluster XC1 
(PL 4) is partially emerging from t he molecular c o re wh ile XC2 is younger and still embed- 
ded. As discussed by REFL08 and iPoulton et al.1 (120081 ). it is not clear that star formation 
here has been triggered by the expanding HII region. 

Since the X-ray sample provides an estimate of the total stellar population of each clus- 
ter, w e can roughly es t imate the star formation efficiency (SFE) of molecular material in the 
RMC. IWilliams et al.l (119951 ) estimate the molecular masses of 70 clumps in the RMC using 
an empirical CO-to-H2 ratio Xrmc — 1 1 x 10 20 H2 molecules cm~ 2 (K km s -1 ) -1 . Mass es- 
timates are made assuming either local thermodynamic equilibrium or gravit ational bindin g 
of gas motions within clumps. If we assume an average stellar mass of 0.4 M (jKroupall2002l ). 
the SFE M star /M star+gas is approximately 15% in region XB and 11% in region XC, adopting 
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the LTE clump masses. If the gravitational clump masses in lWilliams et all (119951 ) are used, 
the SFE becomes 24% and 32% in XB and XC, respectively. These values are similar to the 
SFEs inferred for nearby embedded rich clusters : NGC 2071 (SF E~12%) and NGC 202 4 
(33%) in the Orion Molecular Cloud Lada et alJll99lh . p Oph (9%: Iwilking fc Ladalfl983L 
and Mon R2 (25%; IWolf et alJll99oh . As noted by RL08, the XC (~PL 4) cluster lies at the 
high end of SFEs; this may indicate that on-going star formation has significantly depleted 
gas in the region. 



4.8. Clustered vs. Distributed Star Formation 



Infrared observations have established that 70 — 90% of stars in GMCs form in embed- 



ded clusters (ILada fc Ladal 120031 ) . Unclustered PMS stars have been established mainly in 
the Orion Molecular Cloud complex, but even there it has been difficult to quantify thei r 
population a nd st ar formation processes (IStrom et all Il993l ; lAllenl Il996l ; ICarpenterl 120001 ) . 
Allen et all (120071 ) present new survey results from Spitzer maps of nearby star forming 
clouds, and many show a significant distributed population. Distributed star formation in 
more distant clouds like the RMC is largely unstudied due to the heavy contamination from 
Galactic field stars outside of dense clusters. Modern astrophysical models of star formation 
in GMCs suggests that i ndividual stars may commo nly form asynchronously in supersonically 
turbulent cloud clumps (IMac Low &: Klessenl 12004 ) . 



Because our X-ray identified young stellar population suffers little contamination, we 
are able to examine the fraction of stars that are formed both in clusters and at low stellar 
densities in the RMC. The total number of sources in the three fields is 553, of which 303 lie 
within the contours we use to define clusters XA, XB, and XC ( Table ED . Among the 250 stars 
distributed outside the clusters, at most 10% of the total detections could be Galactic and 
extragalactic contaminants given the short exposure time and the presence of the molecular 
cloud (§ 12.21) . This implies that ~35% of stars could be formed in a distribut ed fashion 
throughout the RMC reg ion, consistent with estimates from infrared studies (e.g. ICarpenter 



2000l ; lLada fc Ladall2003l ). After correcting the sensitivity variation across the fields (§ 14. 3p . 
the fraction of distributed stars is ~45%, which implies that clustered star formation would 
then be the governing mechanism in the RMC but distributed star formation is rather 
important. 

However, further interpretation of this fraction should proceed with caution. There are 
at least three sources of complexity here. First, our definition of cluster boundaries based 
on smoothed stellar density maps (Figure [6]) is arbitrary, which are about 0.8 stars arcmin -2 
for the X-ray-sampled stellar density. In comparison, in REFL08 the cluster boundaries 
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are defined by 0.2 stars arcmin -2 contour levels in the nearest neighbor density map of IR- 
excess sources. They find that 60% stars in the complex are in clusters. The clustering 
fraction reaches 86% if only the stars within the molecular cloud are considered. Besides the 
differences in the field coverage of the X-ray and IR surveys (see also § 15.11) . we could get 
different distributed fractions with different, reasonable choices of cluster boundaries. 

Second, a significant number of lower mass PMS stars (M < 0.5 M ) is missing from 
our survey due to short X-ray exposures (§ 12.31) . While there may be no preference for lower- 
mass stars to be clustered or distributed, it is possible that hi gher-mass stars preferen tially 



form in denser clusters from primordial mass segregation (e.g.. iBonnell fc Davieslll998l ). 



Third, the distributed stars may have been born in a clustered environment but have 
moved into a wider region. Numerical simulations show that the formation process of a 
cluster of stars is highly dynamic a nd chaotic; stellar encounters violently eject some of the 
young stars from multiple systems (IBate et al.l 120031 ) . Sparse clusters in particular can lose 
~ 50% of their stars over a few Myrs. Thus the observed distributed populations may not 
have formed in situ. 



5. Discussion 

5.1. Comparing X-ray and infrared surveys 

An underlying theme to this and similar studies concerns the relationships between 
X-ray and infrared surveys of star forming molecular clouds to uncover young stellar popu- 
lations. Each as advantages and disadvantages. With instruments like FLAMINGOS, NIR 
surveys can cover large cloud areas to faint magnitudes. These catalogs can include virtually 
the entire cloud population at low to moderate extinctions (A v < 20), but often suffer bad 
(often factors > 10) contamination by foreground and background field stars. Used by them- 
selves, PMS samples are restricted to stars with NIR-bright inner disks, usually Class II and 
lightly obscured Class I systems. The resulting NIR-excess samples can constitute only a 
small fraction (< 10%) of the total stellar populations which are often dominated by systems 
with cool or absent disks. With space-based instruments like Spitzer's IRAC, MIR surveys 
also cover large cloud areas with very high sensitivity and penetration to Ay ~ 100. As with 
NIR studies, MIR observations are dominated by field star populations and samples are re- 
stricted to MIR-excess systems. The NIR-excess stars are recovered and new subpopulations 
emerge: embedded Class I and some Class systems, low luminosity PMS brown dwarf sys 



terns, and a few old er systems with transitional MIR-excess disks (e.g.. lBouwman et al.ll2006 



Fur lan et al.l 120071 ). Diskfree Class III members are catalogued but cannot be distinguished 
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from contaminants. 

X-ray surveys have very different characteristics because they trace magnetic reconnec- 
tion flaring rather than stellar and disk bolometric emission. Older Galactic field stars with 
much lower flaring levels are nearly absent, and most extragalactic contaminants are readily 
removed by the absence of stellar counterparts. X-rays can penetrate Ay > 100 but with 
reduced sensitivity. X-ray studies also give estimates of the total stellar populations, and 
measurements of absorption for each star individually. While some Class I and II systems are 
recovered, the main advantage of X-ray observations is that disk-free Class III PMS members 
are readily detected, subject to the sensitivity limitations which correlate with stellar mass. 
Since many clusters are dominated by Class III stars, the X-ray populations can considerably 
exceed the NIR- or MIR-excess populations. A significant drawback is areal coverage: the 
XMM — Newton telescope has a ~ 30' field but suffers confusion in crowded fields, while 
the Chandra X-ray Observatory has a 17' field of view with subarcsecond resolution. This 
limitation can be somewhat mitigated with mosaics as described here. The combination of 
X-ray and infrared studies are more effective at uncovering PMS populations in molecular 
clouds than any one method. We find here that the X-ray samples are quite free from con- 
tamination and can be 10 times richer than the NIR- or MIR-excess samples (Table E]). This 
greatly assists in tracing the spatial structure of star formation in the cloud. Subsequent 
infrared photometry then gives the evolutionary classification, absorption and approximate 
mass of each X-ray-selected star without bias towards luminous disks (Figures [TOlfTll . 



5.2. RMC Young Stellar Clusters 

There is good agreement between our X-ray survey and REFL08 NIR survey (and 
other surveys). Taking this combined X-ray /infrared approach, we emerge with a number of 
results on the PMS population in the portion of the RMC examined with Chandra. First, 
we confirm past findings that star formation is not spatially uniform in the cloud, and we 
estimate that about 1/3 of the stars lie in unclustered environments (§ 14.81) . T his supports a 



large body of infrared-based survey work giving similar results in other clouds (ILada fc Lada 



2003|). 



Second, we estimate the total population in the examined portion of the cloud to be 
about 1700 stars down to the substellar limit (Table El). This estimate is based on an 



extrapolation of the empirical PMS L x -mass correlation (jTelleschi et al.ll2007l ) and knowledge 
of our X-ray sensitivity limits (§ 12.31) . This is similar to the total population of ~ 2000 stars 
found for the central NGC 2244 cluster which illuminates the Rosette Nebula HII region 
(Paper I). 
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Third, in contrast to the single concentrated OB cluster at the center of the Rosette 
Nebula, the RMC population is divided into three smaller clusters which themselves can be 
divided into sparse subclusters. Portions of X-ray clusters XB and XC are associated with 
well-established infrared clusters PL 2 and PL 4 (§ |4.5j) . The X-ray stellar distributions and 
absorptions provide more detail on these structures. XC is richest with a total population 
around 800 stars; subcluster XC2 is more deeply embedded with younger stars than XC1. 
Similarly in cluster XB, group XB2 has more protostars than XB1. These structures are thus 
inhomogeneous and not coeval. We caution that, our current identifications of structural 
components are from Figure [6] (and Figure [7]) rather than Figure H2 as the addition of more 
sources is critically needed for substructure identification. While some substructures were 
cross-identified in the infrared (see Table [5]), the weak substructures certainly require further 
confirmation regarding the details of their locations and sizes. 

Subcluster XB2 = PL 2 is sparse but contains the highest-mass star in our examined 
region of the RMC, the heavily obscured, infrared-excess 09/BO star RMCX #89 = IRAS 
06306+0437 (§ 13.21) . There is only one lower-mass protostar, RMCX #72, lies at a projected 
distance of 0.3 pc to RMCX #89. The absence of a substantial cluster surrounding this 
massive protostar has interesting implications for star formation at the edge of HII regions 
where triggering may be active. The system cannot realistically be dynamically ejected from 
a more crowded region, as it is still embedded in its dense dusty envelope. One possibility is 
that the massive star formed in isolation. We recall that HD 46223, one of the two 04 stars 
in the unobscured NGC 2244 cluster lies in a puzzling location near the southeast boundary 
of the nebula and, based on our Chandra study, has no lower mass companions (Paper I). 
IRAS 06306+0437 may be like HD 46223 seen at an earlier phase. Another possibility is 
that the massive star has formed before (except for RMCX #72) an accompanying cluster 
of l ower mass stars. If thi s scena rio applies, the star formation sequence is opposite to the 



one 



Feigelson fc Townsleyl (120081 ) inferred for the embedded W3 Main clusters, where the 
OB stars appeared to have formed after, not before, an extended cluster of lower-mass PMS 
stars. 

While it is reasonable that the density of obscured X-ray sources decreases significantly 
towards the HII region ionized by the NGC 2244 cluster, it was surprising that the new lightly 
obscured cluster XA appeared in the X-ray survey in a region with no molecular material 
(§ 14.41) . Infrared surveys likely missed this structure due to the paucity of infrared-excess 
stars and confusion by the bright PAH and atomic nebular emission, a stellar enhancement 
associated with XA can be discerned in 2MASS, FLAMINGOS, and Spitzer stellar distri- 
butions. Its location and low fraction of IR-excess sources suggest that it may have formed 
through an early epoch of triggered star formation produced by the expanding Rosette Neb- 
ula HII region. It appears similar to the small NGC 2237 cluster on the western side of the 
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Rosette Nebula (ILil l2005l . REFL08). A detailed Chandra study of NGC 2237 is presented 
in Paper III. 



5.3. Star Formation Modes in the RMC 



While Paper III will present a more comprehensive discussion of star formation in the 
Rosette complex, some immediate inferences can be made concerning the RMC st ructures 



1964 


)• 


Cox et al. 


(1990) 



the compression and triggering proce ss described by lElmegreen fc Ladal (119771 ). However 



neither REFL08, iPoulton et al.l (120081 ) nor the present study support this model in the XB 
and XC region of the RMC. No clear spatial sequence of stellar ages is seen. It is plausible, 
however, that cluster XA was triggered by the expanding Rosette HII region. At present it 
lies in (or, at least, projected upon) the ionized region rather than in the molecular cloud. 
It has the smallest fraction of infrared-excess stars of the clusters examined here, and hence 
was likely the earliest formed. 
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Table 1. Chandra Catalog: Primary Source Properties 
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Note. — Column 1: X-ray source sequence number, sorted by RA. * denotes the sources included in the complete sample of X-ray detections after sensitivity correction (see § |4. 3i ■ 
Column 2: IAU designation. Columns 3,4: Right ascension and declination for epoch J2000.0. Column 5: Estimated la random position error computed as the standard deviation of 
the local PSF inside extraction region divided by the square root of the extracted counts. Column 6: Off-axis angle in arcminutcs Columns 7,8: Estimated net counts extracted in the 
total energy band (0.5 ■ ■ ■ 8 kcV), and the average of its Poisson upper and lower la error (including background contribution). Column 9: Background counts extracted (total band). 
Column 10: Estimated net counts extracted in the hard energy band (2 ■ ■ ■ 8 keV). Column 11: Fraction of the PSF (at 1.497 keV) enclosed within the extraction region. A PSF 
fraction significantly below 90% usually indicates that the source lies in a crowded region. Column 12: Photometric significance computed as the ratio of the net counts to the upper 
error on net counts. Column 13: Log probability that extracted counts (total band) arc solely from background based on Poisson statistics. Some sources have Pg values above the 0.1% 
threshold that defines the catalog because local background estimates can rise during the final extraction iteration. Column 14: Source anomalies: g — fractional time that source was 
on a detector (FRACEXPO from mkarf) is < 0.9 ; e — source on field edge; p — source piled up; s — source on readout streak. Column 15: Variability characterization based on the 
Kolmogorov-Smirnov statistic (total band): a — no evidence for variability (0.05 < Pks)> ^ — possibly variable (0.005 < PkS < 0.05); c — definitely variable (PkS < 0.005). No value is 
reported for sources with fewer than 4 counts or for sources in chip gaps or on held edges. Column 16: Effective exposure time: approximate time the source would have to be observed 
on-axis to obtain the reported number of counts. Column 17: Background-corrected median photon energy (total band). 



Table 2. Chandra Catalog: Tentative Source Properties 





Source 




Position 








Extracted Counts 








Characteristics 






RMCX 


CXOU J 


«J2000 


'5.12000 


Err 


e 


Not 


ANet 


Bkgd 


Net 


PSF 


Signif 


logPfl 


Anom Var 


EffExp 


Me.dE 


# 




(deg) 


(deg) 


(") 


(') 


Full 


Pull 


Full 


Hard 


Frac 








(ks) 


(keV) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


0) 


(10) 


(11) 


(12) 


(13) 


(14) (15) 


(16) 


(17) 


2* 


063232.91+043437.2 


98.137138 


4.577017 


1.4 


8.4 


5.9 


3.4 


2.1 


0.6 


0.90 


1.5 


-2.8 


g... 


13.0 


0.7 


12 


063247.43+043150.2 


98.197650 


4.530632 


1.3 


7.9 


4.5 


3.0 


1.5 


4.0 


0.90 


1.2 


-2.3 


a 


17.6 


2.1 


16 


063249.22+043425.9 


98.205124 


4.573882 


1.1 


7.8 


4.9 


3.4 


3.1 


3.7 


0.90 


1.2 


-1.8 


a 


33.1 


7.8 


26 


063254.45+043724.6 


98.226879 


4.623523 


0.9 


6.3 


4.6 


3.0 


1.4 


3.1 


0.90 


1.3 


-2.6 




35.6 


3.4 


29 


063255.82+043419.1 


98.232618 


4.571979 


1.1 


7.2 


5.2 


3.4 


2.8 


0.0 


0.89 


1.3 


-2.1 


.... a 


34.6 


1.3 


33 


063258.15+043003.7 


98.242303 


4.501051 


1.1 


6.6 


3.3 


2.5 


0.7 


3.5 


0.89 


1.1 


-2.3 


.... a 


15.9 


5.2 


36 


063300.61+043030.7 


98.252556 


4.508536 


0.9 


5.9 


3.5 


2.5 


0.5 


0.0 


0.90 


1.1 


-2.7 


.... a 


17.6 


1.3 


53 


063308.34+044203.4 


98.284762 


4.700962 


1.0 


7.3 


4.6 


3.2 


2.4 


2.4 


0.89 


1.2 


-1.9 


a 


33.9 


1.6 


57 


063308.52+043634.3 


98.285532 


4.609538 


0.3 


4.(1 


4.7 


3.0 


1.3 


0.0 


0.90 


1.3 


-2.7 


g... 


23.3 


0.9 


58 


063309.33+043804.3 


98.288891 


4.634532 


1.1 


6.1 


6.9 


3.7 


3.1 


5.9 


0.89 


1.6 


-2.8 


.... a 


35.4 


3.3 



CO 



Note. — Tablc[^]is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. The column 
descriptions arc the same as in Tabic 1. 



Table 3. X-ray Spectroscopy for Brighter RMCX Sources 





RMCX a 








Spectral Fit b 






X- 


ray Luminosities 




Notes d 


Seq 


CXOU J 


Net 


Signif 


\ogN H 


kT 


log EM 


logL s 


logLh 


logih.c 


log L t 


log L t , c 








Cts 




(cm -2 ) 


(keV) 


(cm" 3 ) 






(ergs s ) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


5 


063239.50+043628.9 


12.3 


2.6 


-0.9 21.8 +0.3 


2.2 


-0.4 53.4 +0.3 


29.67 


29.98 


30.02 


30.15 


30.42 




6 


063239.51+043124.9 


11.7 


2.3 


22.1 


2.0 


53.7 +0.2 


29.72 


30.14 


30.21 


30.28 


30.65 




7 


063240.73+043653.3 


14.4 


2.9 


-0.5 22.2 +0.4 


7.0 


-0.3 53.6 +0.7 


29.48 


30.48 


30.55 


30.52 


30.75 




9 


063245.25+043206.5 


21.4 


3.6 


20.0 


-0.05 0.2 +0.09 


-0.3 53.8 +0.2 


30.45 


27.47 


27.47 


30.45 


30.48 




10 


063246.11+043612.9 


10.8 


2.3 


22.1 


3.5 


53.3 


29.28 


29.98 


30.05 


30.06 


30.34 




14 


063248.65+043404.4 


14.6 


2.7 


20.0 


-0.3 1.2 +0.5 


-0.2 52.9 +0.2 


29.77 


29.10 


29.10 


29.85 


29.86 




17 


063249.80+043641.9 


10.2 


2.2 


-0.7 22.2 +0.4 


3.4 


53.2 


29.07 


29.91 


30.00 


29.97 


30.29 




18 


063251.19+043412.1 


29.9 


4.4 


-0.6 21.9 +0.4 


-0.9 1.7 +2.9 


-0.4 53.7 +0.6 


29.86 


30.02 


30.07 


30.24 


30.60 




21 


063251.73+043621.4 


10.4 


2.3 


-0.3 22.7 +0.2 


2.0 


-0.4 53.8 +0.3 




30.06 


30.32 


30.08 






22 


063252.30+043744.3 


15.1 


2.9 


-0.7 22.1 +0.4 


-2.9 4.0 


-0.4 53.3 +0.6 


29.36 


30.10 


30.16 


30.17 


30.43 





Note. — Table [3] is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and 
content. 

a For convenience columns 1 ■ ■ ■ 4 reproduce the source identification, net counts, and photometric significance data from Table [T] 

b All spectral fits used the "wabs(apec)" model in XSPEC and assumed 0.3Zq abundances. Columns 5 and 6 present the best-fit values for the column density 
and plasma temperature parameters. Column 7 presents the emission measure for the model spectrum, assuming a distance of 1.4 kpc. Quantities in italics were 
frozen in the fit. Uncertainties represent 90% confidence intervals. More significant digits are used for uncertainties < 0.1 in order to avoid large rounding errors; for 
consistency, the same number of significant digits is used for both lower and upper uncertainties. Uncertainties are omitted when XSPEC was unable to compute them 
or when their values were so large that the parameter is effectively unconstrained. Fits lacking uncertainties should be considered to merely be a spline fit to the data 
to obtain rough estimates of luminosities; their spectral parameter values are unreliable. 

c X-ray luminosities are presented in columns 8 ■ ■ ■ 12: s = soft band (0.5 ■ ■ ■ 2 keV); h = hard band (2 ■ ■ ■ 8 keV); t = total band (0.5 ■ ■ ■ 8 keV). 
Absorption-corrected luminosities are subscripted with a c; they are omitted when log./V£r > 22.5 cm -2 since the soft band emission is essentially unmeasurable. 

d 2T means a two-temperature model was used. Well-known counterparts from Table [4] are listed here for the convenience of the reader. 



Table 4. Stellar Counterparts 



X-ray Source Optical/Infrared Counterpart and Photometry 

RMCX CXOU J USNO B1.0 B R I 2MASS FLAMINGOS J II K 2MASS Spitzcr/IRAC 

(mag) (mag) (mag) ID ID (mag) (mag) (mag) PhCcFlg 3.6/^m 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 



1 063232.43+043705.7 ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ yj 

2 063232.91+043437.2 ■ ■ ■ ■ ■ • • ■ ■ • • ■ ■ • • • • ■ ■ • • ■ • • x 

3 063236.58+043611.4 ■ ■ ■ ■ ■ ■ ■ ■ • ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ X 

4 063238.08+043250.1 0945-0095483 ■■• •■■ 14.54 ■•• ••■ ■•• ■•• x 

5 063239.50+043628.9 ■ ■ ■ ■ ■ • • ■ ■ • • ■ ■ • • 063239+043629 14.69 13.45 12.54 y/ 

6 063239.51+043124.9 ■ ■ ■ ■ ■ • • ■ ■ • • ■ ■ • • • • ■ ■ • • ■ • • x 

7 063240.73+043653.3 ■ ■ ■ ■ ■ • • ■ ■ • • ■ ■ • • • • ■ ■ • • ■ • • X 

8 063242.33+043217.9 ■ ■ ■ ■ ■ • • ■ ■ • • ■ ■ • • • • ■ ■ • • ■ • • x 

9 063245.25+043206.5 0945-0095525 18.41 16.59 14.86 06324526+0432074 063245+043207 13.68 12.99 12.79 AAA000 yj 
10 063246.11+043612.9 ■■■ ■■• •■■ ••■ ■•• 063246+043610 14.69 14.29 14.06 yj 



Note. — Tablcl4lwith complete notes is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 

Note. — Columns 1 ■ ■ - 2 reproduce the sequence number and source identification from Table [l] and Tabl e [g] For convenience, [MJD95] — Massey, Johnson, &c Degioia-Eastwood 
(1995), [BC02]=Berghofer & Christian (2002), [TASS4]=TASS Mark IV Photometric Survey of the Northern Sky I Droegc ct al. 2006), [VGK85] — Voroshilov ct al. (1985), [BGD93] =Block, 
Gcballe, & Dyson (1993), [MSX]=the complete MSX6C catalogue in the Galactic Plane. 

RMCX24 = [MJD95] 289 qj 

RMCX25 = [MJD95] 278=[TASS4] 712279=[MSX] G206. 6836-01. 9932; V=12.6 W 



RMCX27 
RMCX32 
RMCX42 
RMOXiS 
RMCX54 
RMCX57 
RMCX58 
RMCX65 
RMCX72 
RMCX89 
RMCX95 
RMCX 119 



[MJD95] 265 

[MJD95] 237=[MSX] G206. 8124-02. 0424 
[MJD95] 230 
[MJD95] 259 
[MJD95] 256 
[MJD95] 260 

[MSX] G206. 7216-01. 9384 

HD 259533=[TASS4] 663214; V=10.6, spectral type GO 
[BC02] 135 

IRAS 06306+0437= [MSXJ G206. 7804-01. 9395 

[TASS4] 712310=[VGK85] NGC 2244 +04 199; V=13.2; spectral type A0V 
= [BC02] 136 



RMCX146 = [MSX] G207.0058-01.9412 



RMCX207 = [BGD93] IRg 4 
RMCX209 = [BGD93] IR g 2 
[BGD93] IRS 5 
[MSX] G207.0302-01.8290 
[MSX] G207.0070-01.8037 
[MSX] G207.1335-01.8414 

[TASS4] 663391=[VGK85] NGC 2244 +04 226; spectral type B9V 
[MSX] G207.1564-01.7804 
[MSX] G207.1792-01.7878 

[TASS4] 1472087= [VGK85] NGC 2244 +04 229; spectral type AO: 
[MSX] G207.1118-01.7013 



RMCX210 
RMCX220 
RMCX241 
RMCX278 
RMCX342 
RMCX352 
RMCX357 
RMCX381 
RMCX387 



I 

CO 

I 



Table 5. X-ray Stellar Clusters and Infrared Counterparts 



X-ray structure IR structure 

RMC RA Deo Extent JHK Spitzer IRAS 

Struc (J2000) (J2000) 



XA 


06:32:38 


+04:46:18 


9' x 6' 






XA1 


06:32:24 


+04:50:30 


1.5' x 3.5' 




06297+0453 


XA2 


06:32:25 


+04:39:00 


3' x 1.5' 






XA3 


06:32:38 


+04:46:18 


6' x 3.5' 






XB 


06:33:15 


+04:35:11 


4' x 8' 






XB1 


06:32:54 


+04:37:00 


3' x 3' 






XB2 


06:33:14 


+04:31:30 


2' x 2' 


C 




XB3 


06:33:15 


+04:35:11 


3.5' x 3.5' 


PL 2 


06306+0437 


xc 


06:34:12 


+04:28:15 


8' x 8' 






XC1 


06:34:12 


+04:28:15 


5' x 6' 


PL 4 E 


06314+0427 


XC2 


06:34:24 


+04:22:00 


4' x 4' 


PL 5, REFL 8 E 


06317+0426 



Note. — JH K clusters are from lPhelps fc Ladal fPL 



(REFL l2008al ). Spitzer clusters are from iPoulton et al 



9971) a nd lRoman-Zuniga et al. 



(120081 ) . 
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Table 6. Star counts in X-ray clusters 



X- 


ray population 






NIR 


pop 




MIR 


pop 


RMC 


N x 


N x ,h 


fx,h 


Not 


Nnir 


Nd,NIR 


fNIR 

J r ] ATTR 
^ ti, J V 1 IX 


fNIR 


Nmir 


fMIR 


Struc 








Total 














1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


XA 


91 


27 


30% 


300 


61 


6 


10% 


2% 


32 


11% 


XA1 


6 


2 


33% 


20 


3 





0% 


0% 


2 


6% 


"X" A 9 


4 


1 


25% 




Q 
O 


n 
u 


U/0 


U/0 


n 
u 


U /o 


XA3 


on 


r 




17/0 


100 


23 


4 


17% 


4% 


11 


11% 


XB 


52 


24 


46% 


200 


37 


8 


22% 


4% 


16 


8% 


XB1 


11 


6 


55% 


40 


5 





0% 


0% 





0% 


XB2 


5 


4 


80% 


20 


3 


1 


33% 


5% 


8 


40% 


XB3 


14 


7 


50% 


20 


11 


4 


36% 


20% 


1 


5% 


XC 


160 


66 


41% 


800 


131 


9 


7% 


1% 


62 


8% 


XC1 


67 


22 


33% 


340 


57 


5 


9% 


1% 


17 


5% 


XC2 


17 


13 


76% 


90 


14 


1 


7% 


1% 


16 


18% 


Distributed 


250 


105 


42% 


1100 















Note. — Column 1: X-ray region denned in Table The Distributed values represent the total 
sources in the X-ray fields minus those in the X-ray regions. 

Column 2: Nx = Number of X-ray sources in the region defined by the contours in Figure Eh and 
b 

Column 3: Nx,h = Number of hard X-ray sources {MedE > 2.0 keV) contained in the region 
Column 4: fxh = Nxu/Nx = Fraction of hard X-ray stars with respect to all X-ray stars. 
Column 5: Ntot = Estimated total stellar population in the region based on the X-ray sample; see 

Column 6: Nnir = Number of NIR stars in the X-ray region from lRoman-Zuhiga et al.l ((2008a) 
Column 7: Nj ntr = N umber of X-band excess disky stars in the X-ray region from 



Roman-Zuhiga et al.l (|2008al ) 
Column 8: f™* R = N d>NIR /N NIR 



Fraction of i^-band excess disky stars with respect to the 
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NIR sample in the X-ray region 

Column 9: f^ f = N^NiR/^tot = Fraction of ET-band excess disky stars with respect to the 
estimated total stellar population in the X-ray region 

Column 10: Nmj r = Number of MIR-band excess disky stars in the X-ray region from 



Poulton et al.1 (j2008h 



Column 11: f^fji = NMin/Ntot = Fraction of MIR-band excess disky stars with respect to the 
estimated total stellar population in the X-ray region 
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Fig. 1. — A 1° x 1° MSX 8.3 /mi image overlaid with 12 CO(J=1-0) emission contours 
(jrleyer et al. 2006 ). outlining the distribution o f the ISM in th e RMC . The small squares 
with sequence numbers denote the locations of iPhefps fc Ladal (119971 ) IR clusters (except 
cluster PL7, which is off the MSX image). The large boxes outline the ACIS-I 17'xl7' 
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Fig. 2. — A Chandra mosaic of the RMC region (0.5-8 keV) overlaid with source extraction 
regions derived by ACIS Extract from appropriate Chandra PSFs. Note that these PSFs 
grow larger in size with increasing off-axis angles. Sources covered by different pointings are 
represented by polygons of different color (red: ObsID 1875; green: ObsID 1876; blue: ObslD 
1877). Note that some sources have two extraction regions as they lie in the overlapping 
regions of two observations. 
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Fig. 3. — The variable X-ray light curve of ACIS source #119. The K - S test finds a 
significant variation (Prs < 0.005), but no flaring characterized by fast rise and slow decay 
as commonly seen in PMS stars is present. 
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Fig. 4. — The X-ray spectrum and the spectral fit for the 07 star ($464). The best fit 
adopts a two-temperature thermal plasma model, with soft kT and a low absorbing column 
N H ~ 4.1 x 10 21 cm" 2 . 
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Fig. 5. — Three-color Spitzer /IRAC+MIPS image of the neighborhood region of source #89. 
Red, green, and blue represents 24/im, 8/im, and 3.6/im emission, respectively. Overlaid 
symbols are the X-ray-selected stars that are Class I (magenta circles), Class II (green 
diamonds), and Class III (black crosses) based on their NIR colors. 




Fig. 6. — (a) The stellar surface density (log^stars arcmin" 2 ) map for all RMC sources 
smoothed with a 3' radius kernel, (b) The same as (a) but using a 2' radius kernel, (c) The 
same as (a) but for the unobscured population, (d) The same as (b) but for the unobscured 
population, (e) The same as (a) but for the obscured population, (f ) The same as (b) but for 
the obscured population. The density scaling is the same between maps of the unobscured 
population and obscured population for a fair comparison. 
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I 




Fig. 7. — Surface density contours based on 10th neighbor distances for the Chandra RMC 
sources. The green contours are densities for all candidates. The blue lines are contours for 
unobscured sources (medE < 2.0keV), and the red lines are contours for obscured sources 
(medE > 2.0keV). All contours follow the levels of density for obscured sources, starting at 
0.057 stars arcmin -2 with steps of 0.45 stars arcmin -2 . 



Fig. 8. — The combined Chandra exposure maps for the RMC observations. The unit is cm 2 
s. Note the enhanced sensitivity in the regions where individual observations overlap. 
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Stellar Surface Density (Log #stars/arcmin 2 ) 




06 h 35™ 00 s 34 m 30 s 00 s 33 m 30 s 00 s 32 m 30 s 

RA (Hour) 

Center: R.A. 06 33 30.44 Dec +04 37 06.1 



Stellar Surface Density (Log #stars/arcmin 2 ) 




06 h 35 m 00 s 34 m 30 s 00 s 33 m 30 s 00 s 32 m 30 s 

R.A. (Hour) 

Center: R.A. 06 33 30.44 Dec +04 37 06.1 



Fig. 9. — The stellar surface density (log^stars arcmin -2 ) map for a sub-sample of RMC 
sources (photon flux > 1 x 10~ 6 count cm~ 2 s _1 ) smoothed with a 3' radius kernel (top 
panel) and a 2' radius kernel (lower panel). The cluster regions and substructure regions 
from Figure O are shown. 
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Fig. 10. — NIR J — H vs. H — K color-color diagram for Chandra stars with high-quality 
JHK photometry (error in both J — H and H — K colors < 0.1 mag) in regions A, B, and 
C as defined in Figure [7J The green circles and red triangles represent sources with and 
without significant i^-band excess (E(H — K) > 2a (H — K)) respectively. The black solid 
and long-dash lines denote the loci of MS stars and giants, respectively, from Bessell & Brett 
(1988). The purple dash dotted line is the locus for classical T Tauri stars from Meyer et 
al. (1997), and the cyan solid line is the locus for HAeBe stars from Lada & Adams (1992). 
The blue dashed lines represent the standard reddening vector with crosses marking every 
Ay = 5 mag. 
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Fig. 11. — NIR J vs. J — H color- magnitude diagram using the same sample and symbols 
as Figure fTUl The purple solid line and dashed line is the 2 Myr isochrone and the 1 Myr 
isochrone for PMS stars from ISiess et al.l (120001 ). respectively. Both isochrones are shown 
here because of the uncertainty in stellar ages. The dash dotted line marks the location 
of Zero Age Main Sequence (ZAMS) stars. The blue dashed lines represent the standard 
reddening vector with asterisks marking every Ay = 5 mag and the corresponding stellar 
masses are marked. 



at ££ s. 



Fig. 12. — Spatial distribution of all the IR counterparts t o our X-ray sourc es. The dense 
molecular cloud is outlined by the 12 CO emission contours (IHeyer et al.ll2006l ). Overlaid red 
polygons are the regions A, B, and C defined by X-ray source densities, and the symbols are 
the X-ray-selected stars classified as Class I (magenta circles), Class II (green diamonds), and 
Class III (black crosses) based on their NIR colors. The 07 star and the Class I protostar 
RMCX ^89 are labeled. The overall distribution of the Class II/I sources is more confined 
to the CO molecular ridge than the Class III stars. 



